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ABSTRACT 

We present an analysis of the 2-150keV spectrum of the transient X-ray pulsar V 0332+53 taken with the 
Rossi X-Ray Timing Explorer (RXTE) in 2004 December. We report on the detection of three cyclotron 
resonance features at 27, 51, and 74keV in the phase-averaged data, corresponding to a polar magnetic field 
of 2.7 x 10 12 G. After 4U 01 15+63, this makes V 0332+53 the second accreting neutron star in which more 
than two cyclotron lines have been detected; this has now also been confirmed by INTEGRAL. Pulse-phase 
spectroscopy reveals remarkably little variability of the cyclotron line through the 4.4 s X-ray pulse. 
Subject headings: pulsars: individual (V 0332+53) — stars: flare — stars: magnetic fields — X-rays: binaries 
— X-rays: stars 



1. INTRODUCTION 

The recurring transient X-ray pulsar V 0332+53 was dis- 
covered in 1983 in Tenma observations made from 1983 
Novem ber to 1984 January (Tanaka Hl983l IMakishima et al. 
1990b). Subsequently, an outburst was r evealed in Vela 5B 
data to have occurred in 1973 summer (ITerrelletalj ri983: 
iTerrell & Priedhorsk\ll984T) . In rapid follow-up observations 
to the Tenma detection by EXOSAT during the 1983 out- 
burst, 4.4 s pulsations were discovered, an accurate position 
was determined, and indications were found for an orbital 
period of 34.25 days (IStella et al.lll985l) . The position and 
subsequent optical/X-ray variability secured the identification 
of the counterpart with the heavily reddened B star BQ Cam 
( Argvl e et alJl!983l). Later inve stigations refined the spectral 
type to 08-9 (Negu eruela et al.ll999l) . 

Analysis of the Tenma data revealed a shape similar to that 
seen in other accreting X-ray pulsars with a flat power law 
and exponential cutoff, a nd also showed evidence of a cy- 
clotron line at -28keV (Ma kishima et"aTlll990bl) . In 1989 
September the source experienced another outburst, this time 
observed by Ginga (Makino 1989). With the energy range 
of the Large Area Counters adjusted to cover the 2-60 keV 
range, cyclotron resonant scatterin g features ( CRSF) were de- 
tected at 28.5 and 53keV ( Makis lfima et alJIl 990ah . At the 
time, this was the fourth accreting X-ray pulsar system to ex- 
hibit cyclotron features. Since then a significant number of 
this class of objects have been shown to have CRSFs, with 
4U 01 15+ 63 setting the recor d for the number of harmonics 
at five (see Heindl et al. 2004, for a review). 
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V 0332+53 went into outburst again in 2004 November , 
reaching a 1.5-12keV intensity of ~1 Crab jRemillarJ2004l) . 
A long serie s of observa tions by the Proportional Counter 
Array (PCA; Uahoda et al.l I2005I) 8 and High Energy X-Ray 
Timing Experiment (HEXTE: lRothschild et al] fl998) aboard 
RXTE were made throughout the outburst. Following up on 
our initial announcement ( Coburn et al. 2005), here we report 
the discovery of a third cyclotron line at ^74keV, along with 
clear evidence for a non-Gaussian shape of the fundamental. 
In §2 we present the observations and the data extraction. The 
phase-averaged and the phase-resolved spectra are analyzed 
and discussed in §3 and §4, respectively. In §5 we summarize 
the results and present our conclusions. 

2. OBSERVATIONS AND DATA EXTRACTION 

Since the luminosity of the source is strongly variable dur- 
ing the outburst, we decided to concentrate on the three 
longest observations during the peak of the outburst as ob- 
served in the soft X-ray band by the RXTE-A11 Sky Monitor 
(ASM; Fig. G). Results fr om the full monitoring data set will 
be presented elsewhere. We combined the PCA and HEXTE 
data taken 2004 December 24-26 (RXTE ObsIDs 90089-1 1- 
05-08G, 90089-22-0 1-00G, and 90089-22-01-01G), giving a 
total of ^58 ks of exposure in the PCA and 14.4ks of dead- 
time-corrected live time in HEXTE cluster B. At the time of 
the selected observations, HEXTE's cluster A was not rock- 
ing due to a single event upset, and we discarded data from 
this cluster 9 . The data were reduced with HEASOFT ver- 
sion 5.3.1, and sp ectral mo deling was performed with XSPEC 
version 11.3.1w ( Arnaudl ll996l) . For modeling the phase- 
averaged spectra (Fig. |2^), PCA data from 3-20 keV and 
HEXTE data from 18-150keV were used, and the HEXTE 
data were rebinned by factors of 3 and 5 between 70-100 and 
100-150keV, respectively. A systematic error of 0.5% was 
assumed for the PCA. 

3. PHASE-AVERAGED SPECTRA 

The X-ray spectrum observed from accreting X-ray pulsars 
is the superposition of X-rays originating in the accretion col- 
umn and the hot spots on the surface of the neutron star. The 
neutron star's X-rays are partially intercepted by the material 

8 See als o|http : / / Iheawww . gsf c . nasa . gov/ docs /xray/xte/pca7| 

9 Cluster A was subsequently rebooted and continues to work normally. 
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FIG. 1.— RXTE ASM 1.5-12keV light curve of the 2004/2005 outburst 
of V 0332+53, binned to 0.5 day resolution. Tick marks denote the times of 
the observations analyzed in this Letter. 



present in the system, such as the stellar wind, giving rise to 
significant absorptio n which ca n exceed 10 23 H- atoms cm" 2 
(see, e.g. JBecker & Worfil2005l for a discussion of the forma- 
tion of the X-ray continuum in the accretion column). 

The 0(1O 12 G) magnetic field close to the magnetic poles of 
the neutron star leads to the formation of CRSFs in the X-ray 
spe ctrum. The obse rver's frame energy of the CRSFs is given 
bv (Mesza rosl 19921) 



£c,„=^(l+yi+2„^sin 2 



1+z ' 



(1) 

(2) 



where n = 1,2,. . . , is the harmonic number, E-p = 11.6keV • 
(B/10 12 G) is the fundamental nonrelativistic cyclotron en- 
ergy (n = 1), 9 is the angle between the wave vector of the 
incoming photon and the magnetic field B, and z ~ 0.3 is the 
gravitational redshift (assuming a 1.4 Mq neutron star with a 
10km radius). The great importance of CRSF observations is 
that they provide the only direct way to measure neutron star 
B-field strengths. 

The lack of adequate theoretical continuum models for ac- 
creting neutron stars necessitates the use of empirical mod - 
els to describe the observations ( Krevkenbohm et al. 1999). 
We use a pow er law modifie d at higher energies by a "Ferm- 
Dirac cutoff" (iTanakal 19861) to fit the PCA and HEXTE data 
simultaneously, including a multiplicative constant as a fit- 
ting parameter to allow for the slight difference (<5%) in the 
flux calibration of both instruments. Photoelectric absorption 
and a strong Fe Ka fluorescence line at 6.35 keV from neu- 
tral iron, with a width of CTFeKa = 0.43 keV and an equivalent 
width of 123 eV are also taken into account. The soft contin- 
uum parameters found (s ee below) are in general agreement 
with earlier results (e.g., linger et alJl992l) . 

While this model fits the overall shape of the continuum 
reasonably well, strong absorption-like features are present 
above 20keV (Fig. |2j?), which we model as lines with Gaus- 
sian optical depth profiles of the form 



t(E) = t c ■ exp I -- 



1 (E-Ec 



rj c 



(3) 



where Eq is the line energy, oc is the line width, and tq is 
the optical depth at the line center. The spectral shape is then 
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FIG. 2. — Cyclotron line modeling of the phase averaged PCA/HEXTE 
spectrum of V 0332+53. (a) Combined PCA/HEXTE spectrum (crosses), 
best fit model (histogram) and unfolded spectrum (dotted and solid lines, right 
y-axis label), illustrating the cumulative effect of the CRSFs on the contin- 
uum (see text for a description of the continuum model). Note that the lines 
are already visible in the raw count spectrum. The inset shows the best-fit 
spectra with one (dotted line) and two (dashed line) line components for the 
fundamental line relative to the best-fit continuum with no lines (solid line). 
Panels (b)-(f) show residuals for models taking different numbers of CRSFs 
into account (see Table 1 for the best fit values), (b) No line; (c) one line at 
25.5 keV; (d) two lines, one at 26.6keV and one at 50keV; (e) two lines, a 
fundamental line that is best modeled with two Gaussian components at 25.2 
and 26.6 keV and its harmonic at 49.9 keV; (f) three lines, the fundamental 
modeled by two components at ~25 and ~27 keV, and two lines at ~51 and 
~74 keV. [See the elctronic edition of the Journal for a coior version of this 
figure.] 



given by mu ltiplying the continuum spectrum by exp(-r(£')) 
jRrev kenboh m et alJ2004 and references therein). 

We start by modeling the strongest feature between 20 and 
30keV with one of these Gaussian components (see Table 1 
for the relevant best fit para meters). This feature corr esponds 
to the CRSF first reported by Maki shima et alJ (fT990a). While 
the 20-30 keV residuals improve significantly compared to 
the pure continuum fit when this component is included, an- 
other strong and remarkably well defined absorption-like fea- 
ture emerges at ~50keV (Fig.|2j;). This feature is the first 
harmonic of the cyclotron line, found at ~ 2 times the en- 
ergy of the fundamental and already tentatively reported by 
MakishimaetaD lQ~990a). Taking these two lines into ac- 
count, weaker r esiduals remain at ~ 25 keV (Fig.[2ji). Similar 
to 4U 01 15+63 dHeindl etalj fl999) and other sources, these 
residuals are caused by the fundamental CRSF having a non- 
trivial shape, which cannot be modeled well with any sim- 
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TABLE 1 

Spectral parameters for the phase-averaged spectrum. 



Parameter 


2c 


2d 


2e 


2f 


N u (10 22 cm- 2 ) 


0.0^-2 




1 5+ 05 

—0 7 


1.6181 


r 


-o.i6;88? 


0.04+0$ 




0.42+°$ 


Ecut (keV) 




+0 7 


13.0^| 


17.2+1-2 


Bfoid (keV) 


6 1+ 00 

-0.0 




7.0+o : | 


8.0+°;? 


Ecyc.iz (keV) 


25 5+ - 

-0.0 


26.3+°;' 


26.6+j]-! 


27 1+ ' 2 


CTcyc.la (keV) 


4 7+0.0 
-0.0 


5 5+ 01 

-0.1 


7.0+0| 


7 6+ - 2 
-0.2 




1 02+ 001 

' -0 1 


i 

1 .Z.*L_Q 


1.38+0- 20 


1 oi +0.08 
-0.06 


£ cy c,lb (keV) 






25.2+°- 2 


25 2+ 2 
—0.2 


o"cydb (keV) 








1 5+°- 3 

—0 5 


7~cyc. lb 






3S+013 
u - JJ -0.04 


34+0 09 
UJ -0.03 


£ cy c,2 (keV) 




50 +0 4 


49 - 9 -a3 


51 3+ ' 4 

J1 - -0.4 


Ccyc.2 (keV) 




6 5+ 05 


6 5+ 05 


g Q+0 6 
°*-0.6 


7~cyc,2 




94+° 09 


1 1Q+0.14 


9 91+0.31 
z - -0.26 


E cy c3 (keV) 








73.71JJ 


Ccyc.3 (keV) 








4 5+ 2 - 5 


'i cyc,3 








t 0+2.6 
J - J -l.l 


xL/dof 


16.24/104 


4.05/101 


1.37/98 


1.04/95 



NOTE. — The columns refer to residuals shown in Fig. 2. Symbols de- 
note the equivalent hydrogen column, TVh, the photon index, T, the cutoff 
and folding energies of the Fermi-Dirac continuum, E cut and £f id, and the 
cyclotron line parameters £ C yc, f C y C , and r cyc (defined in Eq.|3} for fits with 
1, 2, and 3 CRSFs. CRSF parameters labeled la and lb correspond to the 
complex profile of the fundamental line. See text for a description of the 
continuum model. x 2 ed is the reduced \~ and dof is the number of degrees 
of freedom. All uncertainties are 90% confidence levels for one interesting 
parameter (Ax 2 = 2.71). 
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FIG. 3.— Pulse profile and CRSF variability of V0332+53. (a) HEXTE 
(full band, upper profile) and PCA (20-50 keV; lower profile) pulse profile. 
Also shown are the phase evolution of (b) the energy and (c) the depth of 
the fundamental and (d, e) the first harmonic cyclotron line. Error bars are 
at the 90% level for one interesting parameter; dotted lines indicate the best- 
fit values from the two CRSF fit to the phase averaged spectrum. [See the 
elctronic edition of the Journal for a coior version of this figure.] 



pie line shape (see To describe the complex shape we 
add a shallow (r cyCi ib ~ 0.3) and narrow (cr cyc .ib ~ 1.5 keV) 
line component at an energy similar to that of the fundamen- 
tal, iicycib = 25.2^2keV; the shape of this more complex 
line profile is shown in the inset of Fig. [2^. Although low- 
ering the x 2 considerably, the fit is still not acceptable, and 
significant residuals remain (Fig. |2^). Adding another line 
at 73.7+J-3 keV improves the fit significantly to x^ ed = 1.04 
(Fig-Ef, the unfolded spectrum and the consecutive effect of 
the individual CRSFs are shown in Fig.|2t). The F-test prob- 
ability that th is improvement is due to chance is only 3 x 10~ 6 
(however, see Protassov et a 1J2002I concerning the applicabil- 
ity of the F-test in this case). 

To verify the detection of the CRSFs and to confirm their 
independence from the choice of the continuum, we modeled 
the data with other continua commonly used to describe the 
spectra of accreting pulsars. The n egative positive exponential 
continuum (NPEX, Miharal [l995l) . for example, also requires 
the presence of three CRSFs at 27, 51.7, and 74.8 keV, again 
with a non-Gaussian fundamental. Other continua gave simi- 
lar results, as do other prescriptions for the shape of the CRSF. 
The actual choice of the continuum or line shape is thus irrel- 
evant for the detection and characterization of the CRSFs at 
the level of the RXTE energy resolution. We therefore report 
the discovery of the second harmonic of the cyclotron line 
in V 0332+53 with RXT E, confirming our p reliminary initial 
report for this data set (Coburnet al. 2005). This discovery 
has been independe ntly confirmed by sub sequent INTEGRAL 
observations ( Rrevkenbohm et al. 2005). 



4. PHASE-RESOLVED SPECTRA 

The neutron star's magnetic axis is offset from its spin 
axis, not only giving rise to X-ray pulsations but also 
confronting the observer with a complex situation. Due 
to the rotation of the neutron star, the accretion column 
and the hot spots are seen under a constantly chang- 
ing viewing angle. Since the physical conditions are 
likely to be variable within the emission region and since 
the CRSF line shape depends on the viewing angle, we 
expect the X-ray continu um and line parameters to de- 
pend on the pulse phase JArava-G ochez & Harding] 120001: 
lArava & Harding! 1999tllsenberg. Lamb & WangH998l) . This 
has been confirmed by pulse-phase-reso lved spectroscopy in 
a num ber of sources, such as GX 301-2 dKrevkenbohm et all 
|2004[ finding ~20% change i n the fine energy), Vel a X- 
1 ( Rrevkenbohm et al. 2002; La Barbera et al. 2003), or 
4U 0115+63 (Mihara. Makis hima & Nagasel2004l) . 

To obtain phase-resolved spectra, we determined the pulse 
period for each of the three observations using standard tech- 
niques. This approach is necessary, as the accumulated uncer- 
taint y of the only avai lable orbital ephemeris of the neutron 
star JStella et alJl!985l) is too large to allow any meaningful 
orbit correction. The average non-orbit-corrected X-ray pe- 
riod is P = 4.3745 s. After phase-aligning the resulting pulse 
profiles, we determined phase-resolved spectra. Fig.|3ji shows 
the resulting HEXTE and PCA pulse profile of V 0332+53. 

Since the data modes of the observation do not contain a 
PCA data mode with energy resolution below 20keV and 
a sufficiently high time resolution for phase-resolved spec- 
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troscopy, we performed phase resolved spectroscopy using 
the HEXTE data only. Signal-to-noise ratio considerations 
allow a maximum of six phase bins (Fig.|5t). In the resulting 
phase-resolved spectra, the source is detected to ~100keV, 
such that we only model the lower two CRSFs and do not in- 
clude the structure of the fundamental line. The absence of 
soft spectral data forces us to hold T and £ cut at the values 
found from the appropriate phase-averaged spectrum. This is 
justified by the relative constancy of these parameters as seen 
with earlier instruments (Ung er et alJl992|) . Both CRSFs are 
detected with high significance in all pulse phases. Figs.|3]3-e 
show the evolution of the CRSF line energy and depth over the 
X-ray pulse. In contrast to the sources mentioned above, there 
is little to no variation in these parameters, no other spectral 
parameters show any significant variation either. 

5. SUMMARY AND CONCLUSIONS 

We have presented the discovery of a third cyclotron line 
in V 0332+53, making this X-ray pulsar the second, after 
4U 01 15+63, with more than two cyclotron lines, and con- 
firming the identification of the 26keV line as the funda- 
mental CRSF. Similar to sources such as 4U 01 15+63 or 
GX 301-2, the HEXTE spectrum shows that the profile of 
the fundamental line has an asymmetric profile that is shal- 
lower toward lower energies. Such a profile is consistent 
with expectations from the asymmetric shape of the rel- 
ativistic cross section for resonant electron scattering, but 
could also be influenced b y the emission wings pre dicted 
in Monte Carlo simulations ( Arava-Gochez & Harding 2000; 
llsenberg. Lamb & Wang il998> or by a variation of the fi-field 
in the emission region. 

Consistent with the trend seen in most sources with multiple 
CRSFs, the energies of the fundamental and its harmonics do 
not scale like 1:2:3, as predicted by the non-relativistic Eq.|2] 
Rather, the line ratios are slightly smaller, in agreement with 
the prediction of relativistic quantum mechanics (Eq.^. Us- 
ing all three observed line energies (taking E = 26.3 keV for 
the fundamental line) and assuming z = 0.3, for V 0332+53 a 
most probable polar magnetic field of 2.7 x 10 12 G is found. 



Consistency between all three CRSF energies can only be 
reached for 8 > 60°; i.e., photons scattered out of the line 
have an initial direction that is perpendicular to the B-field, 
as is expected from the significant increase of the electron 
scattering cross section with 8 — ► 90°. This result provides 
a much stronger confirmation of the interpretation of the line 
features as CRSFs than is possible in sources where just the 
fundamental is observed. 

Pulse-phase-resolved spectroscopy revealed that 
V 0332+53 shows far less variability in the CRSF pa- 
rameters than many other sources. The origin of the CRSF 
variability in these sources is unclear, although it has been 
pointed out that their complex pulse profiles make it likely 
that different regions of the accretion column or even dif- 
ferent magnetic poles are seen during the X-ray pulse. The 
remarkably simple hard X-ray pulse profile of V 0332+53 
suggests a less complex accretion geometry than in other 
sources, which might also explain the very weak CRSF 
variability. 

In conclusion, our observation of multiple lines strongly 
confirms the interpretation of these features as cyclotron lines. 
For a further interpretation of our results, however, more real- 
istic models for the neutron star continuum and especially for 
the line shape are needed, to include the radiative transfer in 
the line, multiple photon scattering, etc., and thus also predict 
emission wings and strongly asymmetric profiles. Work is 
currently ongoing in our collaboration to compute such mod- 
els for a variety of parameters to make such progress possible 
in the near future. 
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